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The steady state for a system of N particle under the influence of an external field and a 
Gaussian thermostat and colliding with random "virtual" scatterers can be obtained explicitly 
in the limit of small field. We show the sequence of steady state distribution, as varies, 
C , forms a chaotic sequence in the sense that the k particle marginal, in the limit of large N, 

is the /c-fold tensor product of the 1 particle marginal. We also show that the chaoticity 
properties holds in the stronger form of entropic chaoticity. 

O 1 introduction 

00 

\^ I Ever since its introduction by Kac [1] in 1956, the notion of a chaotic sequences has become an 
I important concept in studying many body systems. Chaotic sequences and propagation of chaos 
are the principal tools for passing from many body problems to effective equations. The aim of this 
^ article is to give yet another example of the interplay of chaotic sequences and effective equations. 

In [BDLR] the authors consider a system consisting of N particles moving in a 2 dimensional 
billiard and colliding with convex scatterers that form a dispersing billiard. The particle are subject 
to an external electric field E and a Gaussian thermostat that keeps the kinetic energy fixed. The 



equation of motion between collisions are thus: 



Vi i = 1, . . . , N 



where 

N N 

j(V) = ^^;„ U{V) = J2H' (2) 

i=l i=l 

and Ei is the force exerted on the ith particle by collisions with the fixed scatterers. Very little is 
known about billiards with more than one particle. In particular there is no existence theorem for 
the SRB measure of this system. 

The authors introduced a stochastic version of the above model in which the deterministic 
collisions are replaced by Poisson distributed collision processes. More precisely, in the time interval 



1 



between t and t+dt, the i-th particle as a probability \vi\°'dt of suffering a collision. When a collision 
happens, the velocity of the particle is randomly updated, i.e. if the particle's velocity direction 
before collision is w = |, after collision it will be distributed as K{u' ■ uj)du'. The details of 
the collision kernel K are largely irrelevant. For what follows it will be enough to assume that 
K{x) > for X in a open set W G [— 1, 1]. We note that this stochastic process make sense for any 
dimension d. 

The master equation for this process is given by 

dW{Q,V,t;E) _ ^ giy(Q,V,t;E) ^ T/ ^H^yl.AwiO M f F) 
di " % y'^^^[[E--^v,)W{Q,V,t,E)^ 

N 

+ / K{uj,-uj'^[W{Q,Yl,t;E)-W{Q,V,t;E)]da''~\uj') 

(3) 

where S"^{R) is the m-dimensional sphere with radius R, da^{-) is the uniform surfaces measure 
on S''^{R) and j = J/U as in (2). Moreover 

Q = (gi, . . . ,g7v) , Y = {vi,...,Vi,...,VN) and V[ = {vi, . . . ,v[, . . . ,vn), (4) 

and f ■ = \vi\u^ if Vi = \vi\uJi. Note that the variable Q is not part of the dynamics, i.e., if the 
initial condition iy(Q, V,0) is independent of Q so will be W{Q„V,t). Moreover, if H^(Q, V,0) 
is concentrated on the surface of given energy Uq, that is if 

W{V,0) = 6{U{\) - Uo)F{V,0) 

then so will be the solution of .(3): 

w{y, t- E) = 5{u(y) - Uo)F{v, t- e). 

Finally if -F(V, 0) is a symmetric function so is -F(V, t; E). Thus we will, henceforth, only consider 
symmetric spatially homogeneous solutions concentrated on the surface of energy Uq = N. 

Recall that the A;-particle marginal f^^\vi, . . . ,Vk) of a distribution F/v(V) is defined by the 
equation 



/ y^{v,,...,v,)FN{V)da''^-\V)= [ v{i^,...,Vk)fi;\vu...,Vk)dv,- 



dvk 



where (p{vi, . . . , Vk) is any bounded continuous function on M'"^. Simple computations show that 



y - |Vfc| y5d(iv-fe)-i(|^jv-|Vfc|2) 

where = {vi, . . . , Vk) and V^' = (f^+i, • • • , vn)- A sequence of densities {-F/v}^=i form a chaotic 
sequence with marginal / if for any bounded continuous function ip 



k 

lim / ^(V,)F^(V)rfa'^^-i(V) = / ipiYk) TT f{v^)dv, ■ ■ ■ 



dvk (6) 



It was shown in [2] that for finite time t the master equations (3) propagates chaos, i.e., the 
solution of the master equation (3) forms a chaotic sequence if the initial condition does. More 
precisely, for any bounded continuous function on M*, ip{vi, . . . ,Vk), we have 

k 

lim / ^(Vfe)F(V, t; E)da(V) = [ ^(Vk) TT t; E)dv, ■ ■ ■ dvk 



where 



f{vut;E)= hm fi;\vut;E). 

N^oo 



It was shown in [2] that f{vi,t; E) satisfies the Boltzmann equation 

fivt;E) ^ d \^E_^E.-^t,E))v)f] = {z^ [ K{co ■ u;')[f{v' ,t- E) - f{v,t- E)] da'-\u;') 
at ov I J 

where j{t, E) is given by the self consistent condition 

j{t,E) = - [ vf{v,t;E)dv 
u J 

with 

M = y E)dv. 

The initial condition is given by f{v) = limTv^oo 0). It is easy to see that u is independent 

of time and since we have chosen Uq = N, u = 1. 

Concerning steady states, the situation is far from clear. In [3] and [4] it was shown that a 
steady state -Fss(V; E) exists for the Kac master equation (3) provided that E ^ 0. If = any 
density -F(V) that depends only on the magnitude of the velocities furnishes a stationary state. 
It is, however, an open question whether Fss(V',E) tends to a limiting distribution as -E — > 0. 
Interestingly, assuming the existence of the limiting distribution, it can be computed exactly and 
it is given by 

Fss(V; 0) = i5(f/(V) - iV)F(V) (8) 
Zn 

with 

Fn(V) = i — (9) 

where Zm is the normalization constant 



\2+a 



2+ci 



Zn = / F^(V)da^'^-i(V) . (10) 

For details, the reader should consult [3] and [4]. Thus, the electric field 'selects' the right steady 
state as the electric field tends to zero. 

A similar problem exists on the level of the Boltzmann equation. Again it is possible to show 
that the steady state fss{v', E) for the Boltzmann equation (7) exists and is unique ii E ^ 0. This 
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clearly implies the existence of a steady state current jss{E). In [5] it was shown that, assuming 
that f{v) = \imE-^o fss{v, E) exists and that jss{E) = 0{E), one has 



d 



/(^) = /ile-(v/MH)^+" (11) 



c 

where c and fi are uniquely determined by normalization and the condition u = 1. One easily get 

„ nm , ig^-(i)ir(^) 

which for a = 1 and d = 2 gives 

r f-) 27r / 2 \ 

/X = — )|f ^ 0.65948 c = — r i-)^ 2.83605 



r(i) • 3 V3 

For details the reader may consult [5] where the existence of the small field limit of the steady 
state distribution is proved for d = 1. 

It is now natural to ask whether the distribution (8) is chaotic with marginal / given by (11). 
For the reasons mentioned above, this cannot be deduced from the previous results on propagation 
of chaos since these results do not hold uniform in time. A more serious impediment is the fact 
that the small field limits of the steady states are not known to exist. As explained before, the 
limit as — selects a steady state for the Kac master equation as well as for its Boltzmann 
version. It is far from clear that the selection mechanism is such as to preserve chaoticity. In this 
note we prove that, nevertheless, the distribution (8) is chaotic with marginal (11). 

THEOREM 1.1. Let /Ar(f) be the one particle marginal o/F/v(V) i.e., defined by 



ifiv,) = 1- I F^{Y)da{M') (13) 

and set 

/(^) = (14) 
c 

with the constants given by (12). Then for any bounded continuous function (p{v) 

lim / ip{v)f'^\v)dv = / ip{v)f{v)dv (15) 

and for every k, the k particle marginal fl^\vi, . . . ,Vk) of F{y) satisfies 

k 

lim / if{vi, . . . ,Vk)fN\vi, . . . ,Vk)dvi- ■ -dvk = (p{vi, . . . ,Vk)Y\ fXfJk)dvi ■ ■ ■ dvk (16) 

Af-^oo J^kd Jj^kd 

where, again, (p is any bounded continuous function on 'R^'^. Thus F/v(V) from a chaotic sequence 
with marginal f . 
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2 Proof of Theorem 1.1. 

The following elementary lemma sets the stage for the proof. It will be expressed in terms of the 
probability distribution 

g{w) :-- 



LEMMA 2.1. The following formulas hold for Fn(V): 

N 



r if^) A- |w| 

and 



(V.) : = ^TT^^n / . F^(V,, V^>a^(^-^)"^(V^) (19) 

- iVfci y5d{jv-fe)-i(^;v-|Vfe|2) 

2 + a / N 

■ X 



Proof. Formula (17) follows from (9) and 

A-- = ^J s-'e-^^- , (20) 

r(7) Jo s 

valid for all A > and 7 > by setting 

A= (^H^^"^ , l = dN-l 

and substituting s = t^"*"". The normalization constant Zn given in (10) is then given as 

2 + a 



'dN-l 



f rt'^^-^TTe-(l^^l*)'^"-rfa^'^-i(V) 

r(^l^) Js^'i-HVN) Jo fJi t 



which, using Fubini's theorem and changing variables Vj = yNwj equals 

/ / rre-(l'"^l*^)'^"rfcr^'^-i(W)f 

Jo i^JVd-lfl) f , 



One more variable change r = y/Nt yields 
(2 + «) ""^ " 



^^^^ r f TTe-(l-»^+"rfa^^-i(W)r^^-^- . 

.^+^) Jo Js^'i-Hi)fJi r 
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Hence 

N 



This integral is (18) written in terms of polar coordinates. To see (19) we start with (5) and find 



y - I Vfci r(^^) j5d(iv-fc)-i(^^_|v,|2) Jo * 

Once more, using Fubini's theorem and changing variables vj = \JN — \Vk\'^Wj,i = k + 1, . . . , N 
yields 



2 + a 1 d(jv-fc)-i 



2 X 



Jo fJi >ii ^ ^ ^ 

Changing variables r = \JN — iV^pt yields the expression 





2 + a 1 


1 




iV- 




\'n^)ZN{N-\ 




|2)f 



X 



POO ^ I \ r ^ 

/ ^diN-k)-iTTg / yr ^ ^^^^^ da'^(^-'=)-i(W^)dr . 

which equals 



□ 



LEMMA 2.2. W^e /iawe the following limit 



lim / y^T\g(w,)(m = — 



w/iere 

H ■= \w\'^g{w)dw 

so that 



fdN-l\ ^ (2 + a) 
lim VNT Zn = 

N^OO \2 + a J y/JI 
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Proof. The proof follows with a slight modification from the law of large numbers. We denote 



F{A) := [ T\g{wj)dW , and := / 



— iiYg{w)dw . 



If 



A, = iWe 



|W| 



N 



> e 



then it is a standard estimate that 



From this it follows readily that 



(21) 



(22) 



The set can be written as A^ = Af U Af where 

yl< = {W G R^^ : < /i - and A> = {W G P^'^ • ' ' 



N 



> fi + e} . 



Clearly 



N 



llgiw,)dW < -l=nA 

j=i VP + ^ 



however, to estimate the corresponding integral in the set Af is a little bit trickier, on account of 
the singularity of Fix some number a < \J'^^ independent of so that for sufficiently 
large 

2 



\W\N 



where = [— a, a]'^. Note that the cube C'^'^ = [—a,a]^'^ C Af. By the arithmetic-geometric 
mean inequality 



N 

IwT 



N 



-1/7V 



(23) 



i=i 



so that 



CNd 



N 



N 



N 



l[g{wj)dW 



< 



9[w), \ ^ 1 



It remains to estimate the integral 



L 



N .r 

-Ul,g{w,)dW . 



Here we note that the ball centered at the origin of radius a is a subset of C^'^ and hence 



VA^ 
> - — 



|W| 



and this leads to the estimate 



N 



Collecting these bounds yields 

N 



IW 



< y ^Nd, 



(24) 



and finally 



1 



'N 

W\ 



N 



l[g{w,)dW 



y/jj, — e a 
and using the estimate (21) we find that 



2^ ./jr+^ 



F{A>) (25) 



J^Nd |W| 



N a' 



y/JT^^ e'^N' a e'^N 2^ 



By choosing e = N '^/^ the lemma follows. 



We now turn our attention to . We first need the following elementary Lemma. 



LEMMA 2.3. Let Y = (yi, . . . , y^) he defined by 



\V''\v, 



l<i< k 



yi = Vi k + 1 < i < N . 



Then 



Vi = ya 



N 



0<i< k . 



□ 



(26) 



Moreover the Jacobian determinant is given by 



d{vi, ...,Vk) 



d{yu ■■■.Vk) 



N V Y 



y/c^ /AT^ - Vfc 



N 



2\ 



IV 



k\2 



Proof. We have 



so that 







N- 





iV|Y* 


: 2 




Y|' 


2 



from which (26) follows immediately. To conclude the proof of the Lemma it is enough to observe 
that 

dvi 



N 



Si .; — 



ViVj 



2 I 



□ 
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Let now ^(Vk) be a continuous function on M such that 



sup (piVk) < K . 



Using Lemma 2.3 and 17 we get 



2 + a 



d'YkWgijJi 



i=l 



'.(N-k) 



dY'' --' — -if I 1/1 



ft, 



(27) 



i=l 



Since 



we get 







2 


1 Y'^ 


2 + 




2 



< 1 



\HNiyi, ...,yk)\<K- 



d{N-k)-l 
2+a 



^N~k 



< K' 



for a suitable constant K' . For the last inequality we have used Lemma 2.2. The proof can then 
be completed by applying the Dominated Convergence Theorem if we can show that 



lim HNijji, i/Tv) = (—,...,— ) 



pointwise in W^^. The proof of this fact follows closely that of Lemma 2.2. From (24) it follows 
that 



dY' n 3{yi, 
j=k+i 



N 





Y*: 


bO 


1 Y'^ 


2 + 


Y, 


2 



V I yiTTTT' • • • ' 



2K 



|Y| 



Vk- 



N 



\Y\ 



. yi 



< 



while 



r ^ 

/ dY' n 9iy 



j=k+l 

1 



. N 
'J^'iYfcj 





Y^ 


2 


1 Y*^ 


2 + 


Yk 


2 



lYl 



,yk- 



N 



lYl 



yi 



< 



. sup \ip{\yi, . . . ,\yk) -ip{yu...,yk)\ 

V/^ ~ ^ Ae/]v 

where In is the interval 



N 



N{ii + e) + \Yk\'''\^ N[^^-e) + \Yk\ 



Collecting these estimates we get 
HNiVi,- ■ ■ ,yk) 



. yi Vn 



< 



2 + a 



^ 2+a ■ 



NZ 



N 



. sup \^{\y.u . . . ,\yk) 

_A//i — e AG/jv 





— 1 1 


V/i'' 





IK 



N^, 1 1 



from which Theorem 1.1 follows by choosing e = N s. 

3 Extension and Remarks. 

It is easy to extend the results of the previous section in a couple of interesting directions. We first 
observe that one can give a stronger definition of chaoticity by requiring that given a sequence of 
functions FivfV) on S 



N) , the entropy per particle of this sequence converge to the entropy 



of the 1 particle marginal. More precisely if 



S 



N 



is the entropy of the particles system than 

Sn 



F^(V)logF^(V)cia(V 



lim 



f{v)\ogf{v)dv 



where, as before 



f{v) = lim fii\v). 

N-¥oo 



If this is true we say that the sequence F/v is entropically chaotic, see [6]. 
COROLLARY 3.1. The sequence F]^(y) given by (9) is entropically chaotic and 

d 



lim N-^ 

7V-5>oo 



Fjv(V) logF;v(V)rf(T(V) = log 



/i5 



2 + a 



f{v)\ogf{v)dv (29) 



Proof. We will just report here the minor modification to the proof of Lemma 2.2 needed to prove 
the corollary. We observe that 



xlogx = lim 



^1+5 _ ^ 



5->0 6 



Applying this to (20) we get 



'i< log s'^e-^'ds--i^{-i)A'-^ 



r(7) 7o 

where ip{x) = r'(x)/r(x) is the digamma function. Following the proof of Lemma 2.1, we find 

log^Tv 



dN-1 , /dN-1 



dN-1 (2 + a) 



10 



(30) 



where = c Zn- Using Stirling formula we get that 



N 



dN - 1 fdN - 1 



lim , ^ , 



logc 



d 



2 + a 



Finally we need to compute the integral in the last term of (30). It is very easy to adapt the proof 
of Lemma 2.2 to this expression. Indeed we have that 

, /|w|\ / \og,/jrT-e 

log —= —— < max e, , 

for W G . On the other hand, the function |(logx)/x| is increasing for x < 1 so that from (23) 
we get 



log 



VnJ |w| 



< 



for < 1 



Wi N 



N 



One can now proceed like in Lemma 2.2 with the only difference that a must be chosen such that 

g{w}dw < - 

\w\n ^ 



We finally obtain that 



lim 



/|W|\ ViV-i log^A 



Collecting the above computation we get the first equality in (29). The second equality is imme- 
diate. 

□ 

A last interesting extension regards the first order correction in E. In [4], under the assumption 
that the limit l^'l — )■ exists, it was shown that 

1 / 

F,,(V,E) = ^ Fn{V) + <^^)\v^\HNm + 0{\E\ 



where Vi = IvAui and 



Hn(V) 



c{u) is the unique solution of 



dN-1 



2 + a 



+1 \vi 



1 2+a 



[(Id - /C)c] (w) = u 
where /C is the convolution operator generated by K, that is 



(/Cc) (c.) 



K{uj-uj')c{io')da''''{io'). 
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Since —c{—u) is also a solution if c{uj) is, we have, by uniqueness, that c{u) = —c{—u). As a 
consequence 

5<^-i(l) 

Calling h^^^ the marginal of H^yf, we easily get that 

hPivi, ...,Vk) = -. — IttzVi ■ VvJ'n\vi, ■■■,Vk) 
It is easy to see, from (27), that we can take the limit for — oo on both side and obtain 

k 

lim h^^\vi, ...,Vk) = h{vi) TT f{vk) 

iV— >co -*- 

1=2 

where 

^(^) = J^"^ ■ ^^'/(^) = (2 + «)/^'^/(^) 
Combining the above results we get that the k particle marginal of Fss is 

k \ k 

t^^Hv^ Vk:E)= I l + (2 + a)u-^^ 



lim f^^\v,, ...,Vk;E) = [l + {2 + a)/.^ ■ c{uM Hfivk) (31) 



i=l / i=l 



This is consistent with the results on the Boltzmann equation (7). To solve the steady state 
equation of (7) one as to make an assumption on the form oi jss{E) for small \E. It is natural to 
assume that 

jUE) = tkE + oi\E\) (32) 
where k is the conductivity tensor for the system with 1 particle and energy 1, that is 

c(w) iS)Ujda'^~^{uj). 

5^-1(1) 



Under this assumption one finds that 

Uv, E)=[l + {2 + a)p"-^E ■ c{uj)\v\) f{v) + o{\E\) (33) 

where 

fiv) = ^e-(v^l''l)'^'^ 

with I' and b uniquely determined by normalization and (32). One easily see that the average 
energy of this solution is 



u = \v'^\f{v) dv 



2+ci 
V ^ 2 

SO that, requiring m = 1 we get back the large limit of the one particle marginal of F^^. Clearly 
the first order in E of the /c-fold tensor product of (33) give us back (31). We finally notice that 
the above results tell us that the current per particle at small field for a large system is (2 + «)/i^~ 
times the current of the one particle system, if the energy per particle is 1. 
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